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Abstract 

We consider the following fundamental routing problem. An adversary inputs packets arbitrarily at 
sources, each packet with an arbitrary destination. Traffic is constrained by link capacities and buffer 
sizes, and packets may be dropped at any time. The goal of the routing algorithm is to maximize 
throughput, i.e., route as many packets as possible to their destination. Our main result is an O (logn)- 
competitive deterministic algorithm for an n-node line network (i.e., 1-dimensional grid), requiring only 
that buffers can store at least 5 packets, and that links can deliver at least 5 packets per step. We note 
that O(logn) is the best ratio known, even for randomized algorithms, even when allowed large buffers 
and wide links. The best previous deterministic algorithm for this problem with constant-size buffers 
and constant-capacity links was 0(log 5 n)-competitive. Our algorithm works like admission-control 
algorithms in the sense that if a packet is not dropped immediately upon arrival, then it is “accepted” 
and guaranteed to be delivered. We also show how to extend our algorithm to a polylog-competitive 
algorithm for any constant-dimension grid. 

Keywords. Online Algorithms, Packet Routing, Bounded Buffers, Admission Control, Grid Networks 

1 Introduction 

The core function of any packet-switching network is to route packets from their origins to their destinations, 
but many fundamental questions about packet routing are far from being well understood. In this paper we 
consider one of these questions, namely the competitive throughput network model, introduced by [3], 
Briefly, the model is as follows. The network consists of n nodes (switches) connected by point-to-point 
unidirectional communication links, and we arc given two positive integer parameters, B and c, called the 
buffer size and link capacity, respectively. Executions proceed as follows. Packets are input by an adversary 
over time. Each packet is input at its source node with a given destination node. At each step, each packet 
is either forwarded over an incident link, stored in its current location buffer, or dropped (i.e., removed 
from the system). Storing and forwarding are subject to the constraints that a buffer can store at most B 
packets simultaneously and that a link can cany at most c packets in a time step. These constraints can 
be met for all input sequences since the model allows for packets to be dropped at any time. The routing 
algorithm selects, at each step, which packets are forwarded, which are stored, and which are dropped. The 
goal of the algorithm is to maximize the number of packets delivered at their destination. Since we consider 

‘School of Electrical Engineering, Tel Aviv Univ., Tel-Aviv 69978, Israel. (guy@eng .tau.ac.il). 

1 LIAFA, Universite Paris Diderot, 75205 Paris Cedex 13, France, (moti .medina@liafa. univ-paris-diderot. f r). 

^School of Electrical Engineering, Tel Aviv Univ., Tel-Aviv 69978, Israel, (boazgtau . ac . il). 

§ Supported in part by ISF, MoST and Neptune. 


1 



Ref. 

Dim. 

Comp. Ratio 

Deterministic? 

Range of B, c 

Remarks 

[11, 12, 13] 

1 

0(l°gn) 

/ 

B,c > logn, B/c = n°l 1 l 

immediate from s-t reduction 

[4] 

1 

Oflog 3 n) 

— 

B > 2,c = 1 


[7] 

1 

Oflog 2 n) 

— 

B > 2,c = 1 

FIFO buffers 

[11, 13] 

1 

O(logn) 

— 

B £ [1, logn], c > 1 

also for logn < B/c< n° M 

[12, 13] 

1 

Oflog 5 n) 

/ 

]3,0(logn)] 

preemptive 

[12, 13] 

d 

O(log d+4 n) 

/ 

[3, Oflog n)] 

preemptive 

Theorem 11 

i 

O(logn) 

/ 

[5, 0(logn)] 


Theorem 19 

d 

2°( d ) . log d n 

/ 

[2 d+1 + 1, O(log n)\ 



Table 1: Some results for centralized online algorithms for packet routing. The networks are uni-directional 
grids. In the special case of B = 0 and c > 3, the algorithm in [12, 13] is 0(\og d+2 n)-competitive. 


on-line algorithms, we evaluate algorithms by their competitive ratio, i.e., the minimum ratio, over all finite 
packet input sequences, between the number of packets delivered by the on-line algorithm and the maximum 
number of packets that can be delivered by any (off-line) constraint-respecting schedule. 

It is nearly an embarrassment to find that very little is known about this problem, even in the simplest 
case, where the network topology is the trivial n-node unidirectional line. In this work we provide an 
improved deterministic algorithm for networks whose topology is a d-dimensional grid. 

1.1 Our Results 

Our main result is a centralized deterministic O (log n) -competitive packet routing algorithm for unidirec¬ 
tional lines with n nodes. The algorithm requires buffer size B > 5 and link capacity c > 5. In addition, 
both B and c must be O(logn). We show how to extend the algorithm to d dimensions, where the com¬ 
petitive ratio is 2°^ • log d n, assuming that B,c > 2 d+l + 1. Our algorithm is nonpreemptive, namely, 
packets are dropped only at the time of their arrival (similarly to admission control policies, which “accept” 
or “reject” requests upon arrival). By contrast, preemptive algorithms may drop packets at any time, i.e., 
packets are not guaranteed to reach the destination even after they start traversing the network. The best 
previous deterministic algorithm [12, 13] is preemptive. 

Table 1 provides a summary of our results and a comparison of our algorithm with some previous results 
along various aspects. 
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1.2 Overview of Techniques 

We first explain our approach for the 1-dimensional case. 

The high-level idea is to reduce packet routing in a graph G to circuit switching (or path packing, 
see [15, 6]) in the space-time graph G X T, where T denotes the set of time steps. This so-called space- 
time transformation has been used extensively in this context [5, 2, 7, 16, 11, 12, 13]. To be effective, the 
space-time transformation requires an upper bound on path lengths which does not result in losing too much 
throughput. We use the bound of [12, 13] (which extends [7]), that ensures that the loss is at most some 
constant fraction. After the transformation, we have an instance of online path packing [6, 9]. It is known 
that if the capacities are large enough, i.e., logn < B,c < n°^ l \ then online path packing is solvable 
with logarithmic competitive ratio [6, 11, 12, 13]. We overcome the difficulty that B and c are O(logn) by 
employing a technique called tiling, i.e., partitioning the network nodes into large enough subgrids. Tiling 
has been used in the past [15, 8, 11, 12, 13]; in our algorithm, we use 4 distinct tilings, and work on each 
of them independently. Each tiling induces a new graph called the sketch graph whose nodes are the tiles. 
The capacity of the edges in the space-time graph between adjacent tiles is O(logn) to allow for applying 
0(logn)-competitive path packing algorithms. Path packing algorithms over the sketch graph produce 
sketch paths for accepted packets. Thus, after these preliminary simplifications, we arrive at the sub-task of 
detailed routing, in which coarse sketch paths must be expanded to paths in the original space-time graph. 

Fractional Optimum. Key to our application of the path-packing algorithm is the analysis of Buchbinder 
and Naor [9, 10], which bounds the performance of the algorithm w.r.t. the fractional optimum, which may 
deliver packet fractions. This result allows us to scale buffer sizes and link capacities up and down while 
keeping the competitive ratio under control. 

Combining algorithms. Another central component in the analysis of our algorithm is the combination 
technique introduced by Kleinberg and Tardos [15], Loosely speaking, this technique deals with an admis¬ 
sion control algorithm that is the conjunction of two competitive algorithms, the state of which depends 
only on the requests accepted by both. The technique enables one to prove that the competitive ratio of 
the combined algorithm is the sum (rather than the product) of the competitive ratios of the constituent 
algorithms. 

1.3 Previous Work 

Algorithms for dynamic routing on networks with bounded buffers have been studied extensively both in 
theory and in practice (see, e.g., [1] and references therein). Let us first focus on centralized algorithms 
for d-dimensional grids. We note that while centralized algorithms for packet routing were always relevant 
for switch scheduling, recently the idea of centralization of network functions, including route computation, 
gained substantial additional traction due to the concept of software-defined networks (SDN). See, e.g., [14]. 
The special case of 2-dimensional grids (with or without buffers) is of particular interest as this is the 
underlying topology of crossbars in switches [17]. 

Online Algorithms for Unidirectional Lines. There is a series of papers on uni-directional line networks, 
starting with [3], which introduced the model. In [3], a lower bound of Ul(^/n) was proved for the greedy 
algorithm on unidirectional lines if the buffer size B > 2. Lor the case B = 1 (in a slightly different model), 
an O(n) lower bound for any deterministic algorithm was proved by [7, 4]. Both [7] and [4] developed, 
among other things, online randomized centralized algorithms for uni-directional lines with B > 2. In [4] 
an 0(log 3 n)-competitive randomized centralized algorithm was presented for B >2. In addition, it is 
proved in [4] that nearest-to-go is 0(^/n )-competitive for B > 2. Lor the case B = 1, [4] presented a 
randomized 0(^/n )-competitive distributed algorithm. (This algorithm also applies to rooted trees when 
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all packet are destined at the root.) In [7], an 0(log 2 to)- competitive randomized algorithm was presented 
for the case B > 2. (This algorithm also applies to rings and trees.) In [11], an 0(logn)-competitive, 
nonpreemptive, randomized algorithm was presented. The algorithm in [11] is applicable to a wide range 
of buffer sizes and link capacities, including the case B = c = 1. In [12], an O (log 5 n)-competitive 
deterministic algorithm was presented. The algorithm in [12] is applicable for B,c E [3, log to]. 

Online Algorithms for Unidirectional Grids. Angelov et al. [4] showed that the competitive ratio of 
greedy algorithms in unidirectional 2-dimensional grids is f l(y/n) and that nearest-to-go policy achieves a 
competitive ratio of 0(n 2/3 ’). In [12], an 0(log 6 n)-competitive deterministic algorithm was presented.An 
extension of this algorithm to d-dimensional unidirectional grids, with competitive ratio 0(log d+4 n), is 
presented in [12]. 

For more related results, refer to [13]. 

Organization. The problem is formalized in Section 2. In Section 3 we explain the reduction of packet¬ 
routing to path packing, and the construction of sketch graph. In Section 4 we describe the overall algorithm, 
and in Section 5 we analyze it. Sections 3-5 deal with the 1-dimensional grid (line); extension to the d 
dimensional case is also discussed in Section 6. 

2 Model and Problem Statement 

We consider the standard model of synchronous store-and-forward packet routing networks [3, 4, 7], The 
network is modeled by a directed graph G = (V,E), and by two integer parameter B,c > 0. For the 
most part of this paper, we consider a network whose topology is a directed line of n vertices, i.e., V = 

{n 0 ,.. .,v n - 1 }, E = {(vi-i,Vi | 0 < i < to}. 

Execution proceeds in discrete steps. In step t, an arbitrary set of requests is input to the algorithm. Each 
request represents a packet, and we will use both terms interchangeably. A request is specified by a 3-tuple 
r'i = (a,;. bi,ti), where a, E V is the source node of the packet, b t E V is its destination node, and ti E N is 
the time step in which the request is input. 

In each time step, the routing algorithm removes packets that reached their destination, and decides, for 
each packet currently in the network, including packets input in the current step, whether (i) to drop the 
packet, or (ii) to send it over an incident link, or (iii) to store it in the current node. The selection of the 
action is done subject to the following considerations. 

• If a packet is dropped, it is lost forever. 

• A packet sent from node u over link (it, v) at time t will be located at node v at time t + 1. The link 
capacity constraint asserts that at any step, at most c packets can be sent over each link. 

• A packet stored at node u at time t will be located at node u at time t + 1. The buffer capacity 
constraint asserts that at any step, at most B packets can be stored in each buffer. 

We use the following terminology. A packet n = (a*, b t , tf) is said to be input (or arrive) at a* at time f ; . 
We say that r, is rejected if it is dropped at time t t , otherwise it is accepted. (Our algorithm will guarantee 
that all accepted packets arrive at their destination.) 

Given a set of requests, the throughput of a packet routing algorithm is the number of packets that 
are delivered to their destination. We consider the problem of maximizing the throughput of an online 
centralized deterministic packet-routing algorithm. By online we mean that by time t, the algorithm received 
as input only requests that have been input by time t. By centralized we mean that the algorithm receives 
all requests and controls all packets currently in the system without delay. By nonpreemptive we mean that 
every accepted packet reaches its destination. 
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Competitive Ratio. Let o denote an input sequence. Let ALG denote a packet-routing algorithm. Let 
ALG(a) denote the throughput obtained by ALG on input a. Let OPT(cr) denote the largest possible subset 
of requests in cr that can be delivered without violating the capacity constraints. We say that an online 
deterministic ALG is p-competitive, if for every input sequence cr, |alg(<t)| > - p • |opt(ct)|. Our goal is to 
design an algorithm with the smallest possible competitive ratio. 

3 First Steps 

In this section we present preliminary simplifications we apply to the problem. They include reducing the 
packet routing on a line problem to path packing on grids, and then path packing on sketch graphs. 

3.1 From Packet-Routing on a Line to Path Packing in a Grid 

Let G = ( V., E) denote a directed line with link capacities c and buffer sizes B. The space-time grid 
of G = (V, E) is a directed acyclic infinite graph G st = ( V st ,E st ) with edge capacities c st (e), where 
(i) V st = V x N. Each vertex v G V has infinitely many copies in the space-time grid G st \ namely, vertex 
(v, t) G V st is the copy of v that corresponds to time t. (ii) E st = Eq U E\ where Eq denotes forward 
edges and E\ denotes the store edges. Formally, Eq = {(u,t) — > (v,t + 1) : (it, v) G E , t G N} and 
Ei = {(it, t) -A (it, t + 1) : it G V, t G N}. (iii) The capacity of all edges in Eq is c, and all edges in E\ 
have capacity B. 

The transformation. We transform a request ry = (a,., b t , t t ) for routing a packet in the directed line 
G to a path request rf = ((a,;, t,.). nmfbf) in the grid G st . The coiTectness of the reduction is based 
on a one-to-one correspondence between paths in G st and a routing of a packet in G. Each vertical edge 
(v t , t) -a (vi+i,t + 1) in G st corresponds to forwarding a packet from iy to v l+ \ in step t, and each 
horizontal edge (iy, t ) —> (iy . t + 1) in G st corresponds to storing a packet in iy in step t. 

Embedding in the plane. The naive depiction of G st maps vertex (iy, t) to the point (t. i ) in the plane 
(i.e., the x-axis is the time axis and the y-axis is the “vertex-index” axis). This embedding of G st results with 
a lattice of vertices in which edges are either horizontal or diagonal. We prefer the embedding in which the 
edges are axis parallel, which means that vertex (iy, t) is mapped to the point ft — i , i). In the axis-parallel 
depiction, all the copies of a vertex ty G V still reside in the ith row. However, column j corresponds to a 
traversal of the complete line, stalling at no at time j and ending at iy,_ | at time j + n — 1. 

3.2 From One Grid to Four Sketch Graphs 

Given a grid generated by the transformation above, we apply another transformation to produce a coarsened 
version, called the sketch graph. Specifically, we use tiling. Tiling is a partition of the grid nodes into x i v 
subgrids, where f/, and l v arc parameters to be determined later. We also add dummy nodes to the space- 
time grid G st to complete all tiles. This augmentation has no effect on routing because a dummy vertex does 
not belong to any route between real vertices. 

The tiling is specified by two additional parameters cf> x and <p y called offsets. The offsets determine 
the positions of the corners of the tiles; namely, the left bottom corner of the tiles arc located in the points 

(4> x + i ■ ^h, 4>y + j ' £v), for i, j G N. The algorithm uses four offsets ((f> x , <f y ) G {—ih/ 2,0} x {—i v / 2, 0}. 
We denote these four tilings by T \,..., T 4 . 

Proposition 1. For every vertex (v. t ) of the space-time grid G st , there exists exactly one tiling Tj such that 
(v, t ) is in the south-west quadrant of a tile ofTj. 
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Proposition 1 suggests a partitioning of the requests. 

Definition 2. A request ry = (a,;. h t , tf) is in SW j if the source vertex (a,;, tf) of request r l belongs to the 
south-west quadrant of a tile in the tiling Tj. 

The Sketch Graphs. Each tiling Tj induces a grid, called the sketch graph, each vertex of which corre¬ 
sponds to a tile. The sketch graph induced by Tj is denoted by Sj = ( V(Sj ), E(Sj)), where V(Sf) is the 
set of tiles in Tj. There is a directed edge (si, S 2 ) £ E(Sj) if si f s 2 and E st n (sq x s 2 ) f 0. All edges 
in the sketch graph are assigned unit capacity. 

3.3 Online Packing of Paths 

We use the sketch graphs to solve path packing problems. Intuitively, the path packing model resembles the 
packet routing model, except that there are no buffers, and that each link e may have a different capacity 
c(e). In addition, we generalize the notion of a request to allow for a set of destinations (similar to “anycast”) 
as follows. Usually, the destination of a request consists of a single vertex. If G is a directed graph, then it is 
easy to reduce the case in which the destination is a subset to the case in which the destination is a specific 
vertex. The reduction simply adds a sink node that is connected to every vertex in the destination subset. In 
our setting of space-time grid, the destination subset is a row. Thus it suffices to add a sink node for each 
row (as in [7]). 

Formally, a path request ri in G is a pair (a r , I),) , where a h £ l 7 is the source vertex and I) % C V is 
the destination subset. Let P(ry) denote the set of paths that can be used to serve request r,; namely, every 
path p € P(n) begins in a*, ends in a vertex in I) r , and satisfies some additional constraint (e.g., bounded 
length, bounded number of turns, etc.). Given a sequence R = {ry}j e / of path requests, we call a sequence 
P = {Pi} i£j a partial routing of R if J C / and pi E P (r t ) for every i E •/. The load of an edge e E E 
induced by P is the ratio _ a partial routing of a set of path requests is called a /3-packing if the 

load induced on each edge is at most /3. The throughput of P is simply the number | J\ of paths in P. 

Integral and Fractional Partial Routings. In the integral scenario, a path request is either served by a 
single path or is not served. In fractional routing, a request r, can be (partially) served by a combination 
of paths p\. ■ ■ ■ ■ Pk- Namely, each path p 3 serves a fraction A j of the request, where A j > 0 for all j and 
< L We refer to A j as the flow amount of request r t . The load of an edge e E E induced 

by request r, is the ratio Ylj-eepj A?'/ C ( e )- A fractional solution is /3-packing if the total load on in each 
edge, from all requests, is at most /3. The throughput of a fractional routing is the sum of the flow amounts 
of all requests. Given a fractional routing g, we use \g\ to denote its throughput. Trivially, the maximum 
throughput attainable by a fractional /3-packing is an upper bound on the maximum throughput attainable 
by an integral /3-packing. An optimal-throughput fractional /3-packing can be computed off-line by solving 
a linear program. 

Online Path Packing: Problem and Solution. In the online path packing problem, the input is a sequence 
of path requests R = {r,;},^/-. Upon arrival of a request ry, the algorithm must either allocate a path 
p € P(n) to ry or reject ry. An online path packing algorithm is said to be (a, /3)-competitive if it computes 
a /3-packing whose throughput is at least 1/a times the maximum throughput over all 1-packings. Note that 
for online path packing, we assume that all edges have capacity at least 1 . 

The online path packing algorithm in [ 6 ] (analyzed also by [9]) assigns weights to the edges that are 
exponential in the load of the edges. This load is the load incurred by the paths allocated to the requests 
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req( W—»N) | 

req(S— »N) 

re<7(W—»E) 

req( S->E) 





Figure 1: An a X b grid. The four types of requests (i.e., req(FF— >N),req(S—>N),req(W—,>E),req(S—>E)) 
are depicted by the four arrows. 


that have been accepted so far. The algorithm is based on an oracle that is input r,; and the edge weights, 
and outputs a lightest path pi in P{rf). If the weight of p, is large, then request r t is rejected; otherwise, 
request r t is routed along p,. We refer to the online algorithm for online integral path packing by IPP. The 
competitive ratio of the IPP algorithm is summarized in the following theorem. 

Theorem 3 ([13], following [6, 9]). Consider an online path packing problem on an infinite graph with 
edge capacities such that inf e c(e) > 1. Assume that, for every request ri, the length of every legal path 
in P{ri) is bounded by p max . Then algorithm IPP is (2, log(l + 3 • p majX ))-competitive online integral path 
packing algorithm. Moreover, the throughput of IPP for any request sequence is at least 1/2 the throughput 
of any fractional packing for that sequence. 

Bounded Path Lengths. The load obtained by the IPP algorithm is logarithmic in the maximum path 
length p ma x- This suggests that p max should be polynomial in n. Lemma 4 states that limiting the number 
of store steps per packet by a polynomial in n decreases the fractional throughput only by a constant factor. 

We use the following notation. Given a request sequence R = {r,;} ? ;, let f*(R) denote a maximum 
throughput fractional 1-packing of R, and let /*(i?|p max ) denote a maximum throughput fractional 1- 
packing with respect to R under the constraint that each path is of length at most p max . 

Lemma 4 (after [7]). Let p max = 2n • (1 + f). Then |/*(i?|p max )| > 0.31 • |/*(.R)|. 


3.4 Routing paths across 2-d Grids 

Consider the following special case of routing in grids. Suppose that each path request has a specific source 
vertex which resides on either the south of the west side, and the destination is either the north or the east 
side (i.e., we can route to any vertex on the requested side). For X £ {S, W } and Y £ {N, E), let 
req(4f— >Y) denote the set of path requests whose source is in the X side and whose destination is the Y 
side (see Figure 1). The following claim establishes sufficient and necessary conditions for satisfying such 
path requests. We refer to the routing algorithm used in this case as crossbar routing. 

Proposition 5 (Crossbar routing). The path requests can be packed in the 2-dimensional ax b directed grid 
if and only if \req(W^E)\ + |req(5— >E)\ < a and |req(VF— >N)\ + |req(S— >N)\ < b. 

Proof. The “only if” part is obvious. We now present a distributed algorithm that proves the “if” part. 
Without loss of generality, we may assume that req(5— >N) and i'cc\(\V—tP) are empty. This assumption is 
satisfied by routing such requests along straight paths and giving them precedence over other requests. Thus 
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Figure 2: Satisfying the path requests in the a x b two dimensional directed grid, where a < b. 

we may ignore these lines henceforth, and we are left with the task of routing req (S—yE) and req(VF— >N) 
under the assumption that |req(S'— >E)\ < a and |req(FF— k/V)| < b. 

These requests are served as follows. Order the rows from bottom to top and the columns from left to 
right. Assume, w.l.o.g., that a < b (the case that a > b is solved analogously). 

Requests whose source vertex is in the first a rows or columns turn in the vertex along the diagonal 
emanating from the SW corner. For example, a request in req(FF—»iV) whose source is in row i is routed 
eastward for i hops, and then north for a — i hops (i.e., to the north side of the grid). See Figure 2. 

The requests whose source vertex is in the last b — a columns are routed northward until they reach a 
vertex that does not receive an east-bound path from its west neighbor. Once such a vertex is found, the path 
turns east and continues straight until it reaches the east side of the grid. Indeed, such a right turn is always 
possible because a > |req(iS— kE)| and hence a “vacant row” is always found. □ 

Remark 6 . Proposition 5 extends to the case of capacitated edges assuming all horizontal edges have the 
same capacity and all vertical edges have the same capacity. In this case, the requests can be routed iff the 
number of requests for each destination side is bounded by total capacity of edges crossing that side. 

4 The Packet Routing Algorithm 

We now present the routing algorithm. Pseudo-code is provided in Algorithm 1. The algorithm works as 
follows. First, in lines 1-3, an initial filtering of the requests removes requests if too many requests originate 
in the same space-time vertex (see Definition 9). Then each remaining new request is processed. In lines 
4-6, it is classified as either Near or Far, based on its source-destination distance (see paragraph on packet 
classification). Near requests are routed by the ROUTE-NEAR algorithm, described in Section 4.4. Each 
Far request is associated with the tiling Tj in which its source vertex belongs to a south-west quadrant of a 
tile. Each tiling is processed separately by three procedures: (i) The IPP algorithm, which performs online 
path packing over the sketch graph Sj (line 7). The outcome sketchi is either “REJECT” or a path in a 
sketch graph Sj, i.e., a sequence of tiles from the initial tile to the destination tile, (ii) The INITIAL-ROUTE 
procedure looks for a routing within the SW-quadrant of the first tile of r % : its outcome is either such a path 
denoted initi or “REJECT”. Only IPP and INITIAL-ROUTE may reject a far request. If both procedures are 
successful, then DETAILED-ROUTE is called (line 8). Detailed routing computes a path in the space-time 
graph, i.e., a complete schedule for each packet. In our algorithm, the sketch path for each accepted request 
is computed once and it is fixed, but the future part of a detailed route of a request may change due to the 
insertion of new packets. Therefore, the procedure DETAILED-ROUTE not only computes a path for r, in 
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Algorithm 1 Top-level algorithm for packet routing in the 1-dimensional grid. Code for step t. 


1 

Let R t be a list of new requests, sorted by source-destination distance. 



2 

For each vertex v, let R! t (v) the first B' + d requests in R t whose source is v. 

> filter requests 

3 

for each request r % <G (J ( , R' t (v) do 



4 

if ri £ Near then ROUTE-NEAR(rj) 



5 

else 



6 

Let j € {1,..., 4} be s.t. r* e SW j 


> classify r. 

7 

sketchi <— IPP [Sj,acceptedj,rf) 

> path lengths bounded by p max 

8 

initi -t— INITIAL-ROUTE (acceptedj, rf) 



9 

if sketchi / REJECT and initi f REJECT then 



10 

add r t to acceptedj 



11 

DETAlLED-ROUTE(rj, initi, sketchf) 


> update routes 

12 

else Reject r, 



13 

end if 



14 

end if 



15 

end for 




G st , but may also alter the detailed routes of other requests (without changing the high-level sketch-graph 
routes). 

An important property of IPP and INITIAL-ROUTE is that their state is determined by the requests that 
are actually in the system, i.e., accepted by both. (Rejected requests by either do not affect the state of the 
system.) This property enables us to employ the combination technique of [15]. The listing emphasizes this 
property by explicitly managing the sets of accepted requests for each class (denoted by acceptedj ). These 
sets are arguments of IPP and INITIAL-ROUTE and determine their states. We now proceed to explain the 
algorithm in detail. 

Packet classification. A request n = {a t . 6,. tf) is called near if b- L — a-j < £ v , and far otherwise. We 
denote the sets of near and far requests by Near and Far, respectively. The far requests are further classified 
into four classes denoted by Farj, where Farj = Far n Sw ? . Namely, Farj is the set of far requests whose 
source node is in the SW-quadrant of a tile s in the tiling Tj. 

Link Multiplexing. After classification, there are five classes of requests: one for Near and one for each 
tiling. Routes for each class are computed independently and the final result is the union of routes over all 
classes. This is made possible by partitioning the capacity of each edge in the space-time grid into virtual 
tracks, one track per class. The capacity of each track is | of the capacity of the edge, rounded down, i.e., 
track capacities are B' = \_B/ 5j and d = \cj 5j. (This explains why we require that B,c> 5.) 

Tiling Parameters. Tile side lengths are set so that the trivial greedy routing algorithm is O(logn)- 
competitive for requests that can be satisfied within a tile. Each tile has length and height i v , defined 
as follows. Recall that the maximum path length p inax = 2n • (1 + ^) (cf. Lemma 4). 

Definition 7. We use the following parameters. 

• k = log(l + 3 • p max ) 

• 4 = f01 and l v = I'^1 
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We summarize with the following claim. 

Proposition 8. IfB/c is bounded by a polynomial in n, then the tiling parameters satisfy the following 
properties. 

1. 4 + 4 = O(logrc). 

2. The sum of the edge capacities along each tile side is Of A:). 

3. For each track, the sum of the track capacities along a tile side is at least 6k. 

Proof. Clearly 4 + 4 = Oik). If B/c is polynomial in n, then k = O(logn). The sum of the edge 
capacities along a vertical side is 4 ■ B = Q(k). The sum of the track capacities crossing a vertical side of 
a tile is at least £ v • 5 B' > 6k. The capacities along a horizontal edge is bounded similarly. □ 

Filtering superfluous simultaneous requests with identical sources. Since we do not impose any re¬ 
striction on the requests, it could well be that many requests arrive at the same source vertex in a single time 
step. To deal with that, we use that fact that for each node v and step t, no more than c + B requests can 
leave (v, t ) in any routing. The partition of link capacities for tracks imposes a stricter limitation in the sense 
that within each class, no more than d + If paths can have the same source vertex. 

Definition 9. Given a sequence R of requests, let If denote the subsequence of R defined as follows. For 
each source vertex (a u t, ), choose cl + If packets whose destination is closest to the source node. (If at 
most If + d requests originate at the same node, then all of them are kept in If.) 

Proposition 14 shows that rejecting the requests in R \ R! reduces the fractional optimal throughput only by 
a constant factor. 

4.1 Routing Rules 

The routing at the high level (sketch path) is determined by the IPP algorithm.We now explain the ideas 
behind refining these rough paths (in the sketch graph) into actual paths (in the space-time grid). Throughout 
this section we consider, w.l.o.g., a single tiling Tj. 

Fix a tile s IN t j. We distinguish between the following three types of requests in Farj (we deal with 
Near requests in Section 4.4). 

• Initial requests: requests whose source vertex is in the south-west quadrant of the tile s. 

• Traversing requests: these are requests that enter s from a specific vertex on one side (either west 
or south) and must leave through any vertex of another side (either east or north). The entry vertex 
is determined by a previously-invoked detailed routing, and the exit side is determined by the sketch 
path. 

• Final requests: these are requests whose sketch path ends in tile s. The destination of a final request 
is the north side of s. 1 

Each tile is partitioned into 4 quadrants, denoted NE, SE, SW and NW. We constrain the way requests 
are routed within a tile using the following rules (see Figure 3; no request may cross a thick line). 

1. Initial requests always start in the SW-quadrant and are routed to the north or east side of the SW- 
quadrant along a straight path. The SW-quadrant of each tile is reserved for routing of initial requests. 

2. Traversing requests whose source and destination sides are opposite (e.g., from the south to the north 
side) arc routed along a straight path. 

’The detailed path for n should end in a copy of the destination vertex b, in s. As any path that reaches the north side of s also 
reaches a copy of b;, we use the pessimistic assumption that mw(bi) is the northmost row of tile s. 
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Figure 3: The quadrants of a tile. Traversing requests may not cross the thick border lines. Upper bounds 
on the number of paths in each class that traverse the quadrant sides follow from the path packing algorithm 
IPP. 

3. Paths enter the tile either in the east half of the south side or the north half of the west side. 

4. Paths exit the tile through either the east half of the north side or the north half of the east side. 

4.2 Procedure INITIAL-ROUTE 

Initial routing takes place in the SW-quadrant of the first tile of a Far request. The goal of INITIAL-ROUTE 
is modest: route the request to the boundary of the SW-quadrant. This is done greedily along straight paths 
if possible: it may be the case that some of the edges are already reserved for another request (also routed 
by INITIAL-ROUTE of the same tile). If no straight path is available, INITIAL-ROUTE returns REJECT. 
This means that incoming traffic (of earlier requests) continues uninterrupted along a straight path. Only 
remaining capacity along edges that emanate from a vertex if any, is used for routing the requests that 
originate in (i. j). 

4.3 Procedure detailed-route 

The goal in detailed routing is to compute a detailed path pi in the space-time graph G st given a sketch path 
sketchi in the sketch graph Sj and the initial part of the route initi. The sketch path specifies the sequence 
of tiles to be traversed by the detailed path. In addition, the sketch path specifies the tile sides through which 
the detailed path should enter and exit each tile. Requests that have been assigned a sketch path and an 
initial route must be successfully routed by detailed routing. 

Detailed routing is computed by applying crossbar routing (cf. Proposition 5) to the NW, SE and NE 
quadrants. This routing is computed based on the present requests. As new requests arrive, the future 
portions of the detailed routes may change dynamically so that all requests which are “in progress” will 
reach their destination. Below we argue that crossbar routing indeed succeeds. 

Claim 10. Detailed routing successfully completes the route of each accepted far request. 

Proof By Proposition 5, to ensure successful routing it is sufficient to bound the number of paths that need 
to traverse a quadrant by the capacity of the quadrant side. By Proposition 8, the track capacity of each 
quadrant side is at least 3k. We now prove upper bounds on the number of paths that traverse each quadrant 
side (see Figure 3). The IPP path packing algorithm is a A:-packing over the sketch graph (whose edges have 
unit capacity). It follows that at most k paths traverse each side of the tile. As every request that originates 
in the SW-quadrant of a tile must exit the tile, there are at most 2k paths that traverse each side of the SW- 
quadrant (although their sum is also bounded by 2k). Hence the upper bounds depicted in Figure 3 follow. 
We need to elaborate more on the NE-quadrant because it is also used for routing final requests (i.e., requests 
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that do not exit the tile, but do want to reach its top row). Consider the north side of the NE-quadrant. There 
are at most k traversing requests that wish to exit the tile. In addition, there are at most 2k final requests that 
wish to reach to top row (as each final far request must have entered the tile). Thus, in total there are at most 
3k paths that wish to reach the top side of the NE-quadrant. To summarize, the number of paths that wish 
to reach any quadrant side is bounded by the side’s capacity, and hence by Proposition 5, detailed routing 
succeeds. □ 

Finally, we note that in order for DETAILED-ROUTE to be well defined, we compute it in tiles in column- 
major order, i.e., we start with the bottom tile of the leftmost row and go up, then the bottom tile of the 
second-from left column and go up etc. This ensures that when we reach a tile, all input vertices are fixed. 
We remark that detailed routing can be executed in a local distributed manner; in each time step, each vertex 
needs only to know the initial paths the sketch paths of the incoming packets. 

4.4 Procedure route-near 

Finally, we describe the algorithm for the near requests. The ROUTE-NEAR Algorithm is extremely simple: 
it never stores a packet (i.e., it uses only vertical edges in G st , and gives precedence to older requests). In 
more detail, upon arrival of a request r, 6 Near, the algorithm checks the number of requests already routed 
along the outgoing vertical edge (from ( di,t ) to (a* + l,t + 1)). If this number is less than d, then the 
algorithm routes r, along the vertical path in G st from (a,, t) to (b t . t + (b, — a*)). Note that these edges 
occur in the future, and hence cannot have been saturated by ROUTE-NEAR if the edge outgoing from (a,;, t) 
is not saturated. If there is no free capacity in the outgoing vertical edge, r, is rejected. Note that if r ; is 
accepted, then it guaranteed to reach its destination. 

5 Analysis of Competitive Ratio of the Routing Algorithm 

Our goal is to prove the following theorem for a directed line G of n vertices with buffer sizes B and link 
capacities c, where B,c € [5, log n}. 

Theorem 11. Algorithm 1 is (9 (log n)-competitive with respect to the throughput of a maximum fractional 
routing. 

We translate the problem to a path packing problem over the space-time graph G st . Let ff H , ill) denote 
a maximum throughput fractional routing, and let | ff sL {II) | denote its throughput. Let | ALG {R) \ denote the 
throughput of the online packet algorithm. Theorem 11 follows directly from the following lemma. 

Lemma 12. For every sequence of requests R, \fQ St {R)\ < O(logn) • |ALG(f?)|. 

We outline the proof of Lemma 12. We scale the capacities down by a factor of 0(/r) = 0(logn) in 
the sketch graph. By linearity, this reduces the optimal fractional throughput by the same factor (see Propo¬ 
sition 13). We show that the filtering stage in Line 2 incurs only a constant factor reduction to the optimal 
fractional throughput (see Proposition 14). The filtered requests R' are partitioned into near requests and far 
requests (which are further partitioned into 4 classes, one per tiling). The far and near requests are analyzed 
separately. The analysis of the throughput for far requests builds on the competitive ratio of the IPP algo¬ 
rithm and the INITIAL-ROUTE algorithm (see Claim 15 and Claim 16). By applying the combining analysis 
of Kleinberg and Tardos [15], we show that the competitive ratio for the combined algorithm is the sum of 
the two algorithms (see Claim 17). In Theorem 18, we show that the ROUTE-NEAR algorithm succeeds in 
routing a logarithmic fraction of the filtered near requests. In Section 5.4, the parts of the proof are combined 
together to prove Lemma 12. 
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5.1 Scaling and Filtering 

One advantage of working with fractional routings is that, by linearity, the throughput scales exactly with the 
capacities. Let fg.(R) denote a maximum throughput fractional routing in the sketch graph Sj. Recall that 
the sketch graph has unit capacities. Coalescing of vertices of G s t in each tile results with edge capacities 
that are 0(/c) = 0(logn). Hence, we obtain the following proposition. 

Proposition 13. \f* Sj (R)\ = ^ ■ \fast(R)\- 

Recall that by Definition 9, in the input sequence R! , at most B' + d requests originate in each space-time 
vertex. 

Proposition 14. |/^(i?)|<9-|/^0R')|. 

Proof. Fix a space-time vertex v = (a,, t, ). Let R v (reps. R! v ) denote requests in R (reps. R!) that originate 
in v. Let / = (R). Consider the flow //9. For every vertex v, the amount of flow that originates in v is 

bounded by (B + c)/9 < B' + d. Divert flow from in //9 from R v \ R' v to R' v along shorter paths, to obtain 
a flow g with respect to R' such that \g\ = |/|/9. Since |/^, st (f?')| > \g\, the proposition follows. □ 

5.2 Far Requests 

Two algorithms determine whether a far request is rejected: (i) the IPP path packing algorithm over the 
sketch graph, and (ii) the INITIAL-ROUTE algorithm that deals with routing in the initial SW-quadrant of 
the source tile. We begin by showing that, if invoked separately, each of these algorithms accepts at least a 
constant fraction of the maximum fractional throughout over the sketch graph. 

Let R'j denote the subsequence of requests in R' that are in the class Farj. Suppose we invoke the IPP 
algorithm in isolation over the sketch graph Sj with the input sequence R'-. By isolation we mean that the 
accepted requests are determined solely by IPP. Let 1 1PP^. (/?,'-)| denote the number of requests that are 
accepted by this invocation. 

Claim 15. |/| j ( j R')|< 0 2 t .|ipp S j (^)|. 

Proof By Lemma 4, the restriction of the path lengths by p max only reduces the fractional throughput by 
a factor less than 0.31. By Theorem 3, the IPP algorithm is (2, A:)-competitive, and hence its throughput is 
half the optimal fractional throughput with bounded path lengths. □ 

Let |lNlTlAL-ROUTE(i?' )| denote the number of requests that are accepted by INITIAL-ROUTE if in¬ 
voked in isolation with the input sequence . 

Claim 16. |/^ (R'.)| < 2 • |initial-route(R'-)|. 

Proof A far request must exit the tile in which it begins. The edge capacities in the sketch graph are unit. 
Hence, the amount of flow in fg.(Rj) that originates in each tile is at most 2. On the other hand, if a 
positive amount of flow originates in a tile s, then at least one request starts in the SW-quadrant of s. Hence 
INITIAL-ROUTE(i?' ) accepts at least one request that begins in s. P 

A naive analysis of the requests accepted by the conjunction of the IPP and INITIAL-ROUTE algorithms 
implies that the accepted requests are in the intersection, which might be empty. However, in our algorithm 
the subsequence of accepted requests is determined by both algorithms, and this set of accepted requests 
determines the state of both algorithms. Hence, by applying the combining analysis of Kleinberg and Tar- 
dos [15], the combined competitive ratio is shown to be the sum of the isolated competitive ratios. 
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Claim 17. f* s .(R') < (2 + ^ ■ \alg(R')\. 

Proof. To simplify notation, let A\ denote the IPP algorithm in isolation, A 2 denote the INITIAL-ROUTE 
algorithm in isolation, and ALG denote the combined algorithm. Let I denote the input sequence R!-. 

Consider the execution of ALG with the input sequence I. Let X denote the subsequence of requests 
accepted by the combined algorithm A on input I (i.e., X = A(I)). Let Xj denote the subsequence of 
requests accepted by algorithm Aj during this same execution. Note that X = X\ n X2. Let us rewrite 
fs (I) as a function /*:/—>■ [0,1] in which f*(i) equals the amount of flow assigned to request i £ / 
by fg.(I). We abuse notation and view X also as its characteristic function (i.e., X(i) = 1 if and only if 
i G X). 

Consider the packet routing version in which requests have demands in [0,1]. A demand of 1 corre¬ 
sponds to the situation till now. A demand of zero means that the request does not appear (the input skips 
over this request). A fractional demand means that (i) at most this fraction can be routed by the fractional 
routing, and (ii) the request only occupies this fraction of the capacity of an edge. Fractional demand func¬ 
tions can be added and multiplied. If both functions are integral, then addition corresponds to the union, and 
multiplication corresponds to the intersection. 

Consider the demand function Zg = X + f* ■ (1 — Ag), for i £ {1.2}. First, note that it attains values 
in [0,1] because (i) an accepted request i 6 X satisfies Ag{i) = 1, and (ii) a rejected request i ^ X satisfies 
X(i) = 0. Note also that f* ■ (X + 1 — Ag) is a feasible fractional flow with respect to the demand Zg, 
hence 


i 4 (^)i>iir-(x+i-A,)i| 1 . a) 

We claim that the isolated algorithm Ag on input Zg accepts exactly X. To prove this consider the 
sequence of states of Ag(Zg) and A(X). Let S\ denote the prefix of accepted requests by Ag(Zg) till but 
not including the arrival of request i. Similarly, let X 1 denote the prefix of accepted requests by A{X) till 
but not including the arrival of request i. If S\ = X 1 , for every i, then Ag(Zg) = X. Now we prove that 
S}. = X 1 by induction on i. Before the first request, both sets are empty. The induction step for i G I is 
easy because Ag accepts i when the state is X 1 . On the other hand, if i £ X and Zg(i) > 0, then Ag(i) = 0. 

Let CRg denote the competitive ratio of Ag. Then 

\r s .(Zg)\ < CRg ■ Ag(Zg) = CRg ■ \X\ ( 2 ) 

By Equations 1 and 2 

\\f* .(X + l-Agfh <CRg-\X\. ( 3 ) 

Observe that 1 < {X + 1 — A\) + (X + 1 — Af). Indeed, if i E X, then Ai(i) = ^ 2 ( 1 ) = 1, so both sides 
equal 1. If i ^ X, then A 1 (i) or .4,(2) equals zero (perhaps both), and hence the right hand side is at least 
1. Hence, 

unit < 11/* • (x +1 - A^ih +11 r -(x+i- a 2 )\\i 

< CRg ■ \X\ + CR 2 • |.Xj, 


and the claim follows. 


□ 
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5.3 Near Requests 

In this section we analyze the competitive ratio of the ROUTE-NEAR algorithm with respect to near requests. 
Recall that: (1) A request is a near request if the distance from the source to the destination is at most i v . 
Note that l v = 0(^fp) and B < log(l + 3p max ) = O(logn). (2) The incoming requests are filtered so 
that at most B' + d requests originate in every space-time vertex. 

The following theorem states that ROUTE-NEAR succeeds in routing at least a logarithmic fraction of 
the filtered near requests. This theorem implies that the throughput is at least a logarithmic fraction of the 
optimal fractional routing of the filtered near requests. 

Theorem 18. |ALG(Near)| > • |Near|. 

Proof. It suffices to prove that \ALG(Near)\ > Q • \Near \ ALG(Near)\. Consider the following 

bipartite conflict graph. Nodes on side L are the requests of ALG (Near), and nodes on side R are the requests 
of Near \ M.G(Near). There is an edge (rj, rj ) € L x R if rj is rejected by the ROUTE-NEAR Algorithm 
and the vertical route of rj traverses the source vertex ( aj,tj) of ry. A request r, G L conflicts with at most 
B' + c' requests in each vertex. Hence, the degree of r r in the conflict graph is at most (B 1 + d) ■ i v . On 
the other hand, the degree of rj € R equals d (where d is the capacity of the track reserved for the near 
requests). By counting edges on each side we conclude that 


Y2 deg(r-i) = deg(rj). 

Vi£L rjG-R 


Hence, 


(B' + d) ■ t v ■ \L\ >d-\R\. 


We conclude that 


L 


> 


(B> + d)£ r 




As (B +c )i v _ = O(logn), and the theorem follows. 


□ 


5.4 Putting Things Together 

In this section we prove Lemma 12. We partition the input sequence R' into Near and R!- ; , for j G {1,2,3,4} 
(recall that R'!■ = R' fl Farj). By subadditivity, 

4 

|/£4#)l < \fa*t(Near)] + ^ |/^(^)l, 

3 = 1 
4 

|ALG(i? / )| = \ALG(Near)\ + ^ |ALG(i?' )|. 

3 =1 

In order to bound the ratio \f^, st (R')\/\ALG(R')\, it suffices to separately bound the ratios of the terms. 
Indeed, by Theorem 18 

\fQ S i(Near)\ < 0(logn) ■ \AhG(Near)\. 


(4) 

(5) 
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By Proposition 13 and Claim 17, 


|/^ t (i?')|<0(logn)-|ALG(i?')|. 

Finally, by Proposition 14, |/^ st (i2)| < 9 • \f^ at (R')\. Since ALG(-R) = ALG(i? / ), the lemma follows. 

6 Extension to d-Dimensional Grids 

The following theorem is proved by extending Algorithm 1 for a line network to (/-dimensional grid. 

Theorem 19. For B,c € [2 rf+1 + l, log n], there is a deterministic 2°^ •log''* n-competitive online algorithm 
for the throughput maximization problem. 

(sketch). As in the one-dimensional case, perform a space-time transformation on the (/-dimensional n- node 
grid G to obtain the (d + 1)-dimensional space-time grid G st . Partition G st to 4 x ... x (-d+x subgrids 
(or subcubes). The side length of a subgrid equals l v for directions that correspond to forward steps and 
(/, in the direction that corresponds to store steps. There are two offsets per dimension, resulting with 2 d+l 
tilings. The number of tracks equals the number of offsets plus one (the extra one is for the near requests), 
hence we require that B,c> 2 d+1 + 1. Similarly to the 1-dimensional case, a request is classified as a near 
request if the distance from the source to the destination is at most d ■ l v . Detailed routing within a tile is 
successful by the following observation. Every time a packet cannot turn to the direction that is dictated by 
its sketch path, there is a packet that did turn to its desired direction. Since the number of path emanating 
from each tile is bounded by the quadrant-side capacity, we conclude that every packet will eventually turn, 
if needed, within its quadrant, thus respecting its sketch path. 

Since the link capacity to track capacity ratio is 0(2 d ), this scaling of capacities incurs an 0(2 d ) factor to 
the competitive ratio. The sketch graph is obtained in the same way, with the exception that edge capacities 
are set to ^j-j- (so that that the number of paths that IPP routes out of a (d + 1 (-dimensional tile is at most 
0(log n)). The ratio of edge capacities in G st between adjacent faces of tiles and the capacity of the edge in 
the sketch graph is 0(log n) d . This incurs an additional factor of 0(log n) d for routing far requests due to 
capacity scaling. The routing of near requests succeeds in routing at least a fraction of d ■ log n of the near 
requests. We conclude that the competitive ratio is determined by the fasr requests, and hence the theorem 
follows. □ 

7 Conclusion 

In this paper we presented an online deterministic packet routing algorithm. For the one dimensional grid 
(with constant-size buffers and constant-capacity links), this algorithm closes the gap with the best through¬ 
put achieved by a randomized algorithm. This closes a problem which was open for more than a decade, 
but still leaves open quite a few problems. The most urgent one is to reduce the gap between the upper 
and lower bounds on the competitive ratio. Currently the best upper bound is O(logn) for the line, and we 
are not aware of any no non-trivial lower bound. We note that reducing the upper bound to o(log n) seems 
to require new techniques, as the reduction to online path packing introduces a logarithmic factor in the 
competitive ratio. 

Another important question is to come up with reasonable distributed algorithms. Even though, as 
mentioned above, the SDN model shifts many network operation tasks to the centralized setting, it is very 
interesting to find out what can be done without a central coordinator. 
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A v l°^ //-competitiveness of initial routing 


Lemma 20. Fix a tile s and let q denote its .S’IV quadrant. Suppose that the sources ofm path requests are 
in q. Then path requests are served by the initial routing in q. 

Proof. We restrict attention to rows of q which contain at least B' sources and columns of q which contain 
at least d sources of requests which cannot be routed horizontally. Since all other packets arc trivially routed 
by the algorithm, we may assume w.l.o.g. that there are no other packets with sources in q. 

Let y denote the number of rows that contain a source vertex of an initial request, and let x denote the 
number of columns that contain a source vertex of an initial request that is not routed horizontally. Clearly, 
m < xy(B' + d). On the other hand, detailed routing in q serves yB' + xd requests. We now prove that 
yB' + xd = Fl{y/xy{B' + d)). 

Without loss of generality, assume that d > B'. Thus it suffices to prove that 

/f .g + f-.e = 


We proceed with case analysis. If y < x, then ■ B' + . If ■ d > d, as required. Otherwise y > x. 


We further distinguish between two cases: 

1. If y/x > d/B' , then ^ • B' + • d > • B' > \fd, as required. 

2. If x/y > B'/d, then • B' + ^/| • d > \p^r • d > sfd, as required. 


□ 


Note that, if a single requested is input to a SW-quadrant, then intial routing accepts it. 
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